Geometry versus Entanglement in Resonating Valence Bond Liquids 
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We investigate the behavior of bipartite as well as genuine multipartite entanglement of a res- 
onating valence bond state on a ladder. We show that the system possesses significant amounts 
of bipartite entanglement in the steps of the ladder while no substantial bipartite entanglement is 
present in the rails. Genuine multipartite entanglement present in the system is negligible. The 
results are in stark contrast with the entanglement properties of the same state on isotropic lattices 
in two and higher dimensions, indicating that the geometry of the lattice can have important im- 
plications on the quality of quantum information and other tasks that can be performed by using 
multiparty states on that lattice. 



I. INTRODUCTION 

Resonating valence bond (RVB) states were introduced 
in 1938 by Pauling in organic, and later in metals and in- 
termetallic, compounds [l|. It was extensively studied in 
many body physics ever since Anderson, in 1973, pre- 
sented the idea of using such states to explain the behav- 
ior of Mott insulators [2| . The importance of RVB states 
grew immensely after the proposition of relatingsuch 
states with high temperature superconductivity [1, 
Investigations on the use of RVB states in quantum infor- 
mation has been initiated in recent years. The possibility 
of topological quantum computation by using such states 
have been proposed 5] , and the entanglement properties 
of RVB states in isotropic two- and three-dimensional 
lattices have also been explored @] . 

In this paper, we consider the RVB state on a ladder, 
a pseudo two-dimensional system (See Fig. 1). Ladder 
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FIG. 1: (Color online.) RVB liquid on a ladder. Blue and Red 
balls belong to two different sublattices A and B respectively. 
One particular dimer covering is shown. Singlets are always 
from sublattice A to sublattice B. 

states (obtained by assembling spin chains one next to 
the other to form ladders of increasing width) have re- 
cently generated a lot of interest in theoretical and ex- 
perimental condensed matter physics. Ladder states have 
been proposed in the study of superconductivity in sys- 
tems like doped (VO)2P207 [3] and the study of transi- 
tions from 1-D chains to 2-D spin systems 8] . It has been 
found that ground states of antiferromagnetic Heisenberg 
models (with suitable couplings) on "even ladders" (con- 
sisting of two or more even number of chains) can be RVB 



states @-[T3l- Quantum information concepts like entan- 
glement, fidelity, of different ladder systems have been 
studied (111] , and protocols for high fidelity transmission 
of quantum states through such systems have also been 
presented [T^ . 



Main results 

We study the entanglement properties of RVB states 
on a pseudo-two dimensional ladder Precisely, we con- 
sider ladders of two chains (2 x M ladders, N=2M spins), 
of different lengths. We find that the entanglement char- 
acteristics of an RVB state on a ladder are significantly 
different from that on an isotropic two or higher dimen- 
sional lattice. Bipartite entanglement between neighbor- 
ing sites are of two different types for the case of a ladder. 
While the bipartite entanglement of the neighboring sites 
on the rails of the ladder is insignificant, that of the steps 
(or rungs) is substantial. This automatically suggests 
that genuine multiparty entanglement of the ladder RVB 
is negligible. This is in sharp contrast with the situation 
on an isotropic 2D or 3D lattice, where the single type 
of nearest neighbor bipartite entanglement is negligible, 
while the state is genuine multiparty entangled [6|. We 
reach our conclusions by numerical simulations as well 
as via analytical bounds on the entanglement measures. 
The symmetry of the system leads us to use the singlet 
fraction as a measure of bipartite entanglement, while in 
the multiparty case, we use a recently proposed measure 
called "generalized geometric measure" [ij, [l3| • 

This change in the entanglement properties can be at- 
tributed to the change in the geometry of the RVB sys- 
tem: a ladder is not isotropic in all directions. Our results 
show that this change in geometry has a marked effect on 
the entanglement properties of the system, and hence on 
quantum information tasks possible. The results are ob- 
tained by using periodic boundary conditions. However, 
we have also found the similar trend of bipartite entangle- 
ment, when there is open periodic boundary condition. 
Also the calculations were done for the, physically inter- 
esting, so-called RVB liquid, which we now define. 
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II. RVB LIQUID 

The pseudo two-dimensional lattice under considera- 
tion can be divided into sub-lattices (see Fig. 1), A and 
B, in such a way that all the nearest neighbor (NN) sites 
of any site on sub- lattice A belong to sub- lattice B, and 
vice-versa. Such a lattice is called a bipartite lattice. 
Each lattice site on such a bipartite ladder is occupied 
by a qubit, which can for example be the spin state of 
an electron. For such a system, the RVB state is defined 
as the equal superposition of all the possible dimer (sin- 
glet) coverings that can be formed by nearest-neighbor 
directed dimers. Between site i in sublattice A and site 
j in sublattice B, a directed dimer is defined as 
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where | f) and | l) are respectively the spin-up and spin- 
down states in the 2;-direction, at the corresponding lat- 
tice site. The (unnormalized) RVB state can then be 
written as [ill 



|(ai, 6i)(a2, &2)...(aAr, bpf)) 



(1) 



where and bj are site positions and the summation is 
over all dimer coverings such that the dimers in every 
covering satisfy ai € A and bj € B. Since the dimer 
coverings that constitute the RVB state are formed by 
using only NN dimers, we call the state as an RVB 
liquid. 



III. BIPARTITE ENTANGLEMENT OF RVB 
LIQUID 

We are now ready to study the behavior of nearest- 
neighbor bipartite states for an RVB liquid on a ladder. 
To obtain the bipartite density matrix, pi2, between any 
site (say, 1) in A and one of its nearest neighbors (say, 2) 
in B, we take the partial trace of the whole RVB liquid 
over all sites except 1 and 2: 



Pi2 = TrY2 IV-) (V'l 



(2) 



where Tr^g represents the partial trace over all sites other 
than 1 and 2. Rotational invariance of 1-0), given in Eq. 
(H)), implies that pi2 (and all other reduced density ma- 
trices of I "(/;)) is also rotationally invariant. The only ro- 
tationally invariant two-qubit states are the singlet and 
the maximally mixed state. Therefore, pi2 is the Werner 
state [16] given by 



Pi2{p) ^p\{i,j)) {{i,j)\ 
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(3) 



where the "Werner parameter" p satisfies —1/3 <p< 1, 
and I4 is the identity operator of the four-dimensional 
complex Hilbert space. Bipartite entanglement mea- 
sures, like concurrence 53], are monotonic functions of p. 



Hence, instead of calculating entanglement measures ex- 
plicitly, we will investigate the behavior of p with respect 
to increase of system size. Note here that the Werner 
state is entangled for p > 1/3. 



A. Properties of nearest-neighbor bipartite states 
for periodic boundary conditions 

In this subsection, we will investigate the bipartite en- 
tanglement of two neighboring sites of the ladder with 
periodic boundary conditions. 



1. Bipartite Entanglement: Steps versus Rails 

As mentioned earlier, the nearest-neighbor bipartite 
state of the RVB liquid is a Werner state. By symmetry, 
there are only two types of such bipartite states: along 
any one of the rails, and on a step. We denote the Werner 
parameter p of a Werner state along a rail as Pr, while 
that on a step as Ps (see Fig. 1). 
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FIG. 2: (Color online.) Bipartite entanglement on the rails. 
The decrease of the entanglement parameter pr with increas- 
ing is clearly seen. 



As shown in Fig. 2, the values of Pr consistently de- 
crease with the increasing system size N, and hence bi- 
partite entanglement on the rails decreases with the in- 
creasing TV, similar to the case of an isotropic lattice. In- 
terestingly, however, the BE of the steps increases with 
respect to N (see Fig. 3). 

Therefore, already at the level of bipartite entangle- 
ment, one obtains a trade-off between the two different 
forms of nearest- neighbor entanglement, in the case of 
an RVB liquid on a ladder. Contrast this with the case 
of a square lattice or any other isotropic lattice, where 
the corresponding RVB liquids have only a single type of 
nearest-neighbor entanglement, and its value, quite gen- 
erally, decreases with increasing system size 16|. The nu- 
merical evidence presented above for this complementary 
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The average fidelity of tlie transfer is given by 
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FIG. 3: (Color online.) Bipartite entanglement on the steps. 
In sharp contrast to the isotropic case, bipartite entanglement 
on the steps, as parametrized by ps, increases with increasing 
system size A^. 



behavior of the two types of BE, present in the ladder, 
will later in the paper also be justified by using analytic 
bounds. 



2. Regional Entanglement 

To study the behavior of BE between nearest neigh- 
bors, we introduce a quantity, called "regional entangle- 
ment" (denoted as Pavg) for each lattice site. For every 
lattice site, there are three NN sites, which form three 
nearest-neighbor bipartite entanglements with that site. 
(Note that we are still in the regime of periodic boundary 
conditions.) Regional entanglement at a particular site 
is defined as the average of these three nearest-neighbor 
entanglements, as quantified by their Werner parame- 
ters. Due to periodic boundary conditions, the regional 
entanglement at each site is the same, so that it is a 
characteristic parameter for a given ladder lattice size. 
It is a measure of the regional distribution of bipartite 
entanglement at any lattice site. 

The physical significance of Pavg can be understood by 
using the concept of fidelity of teleportation [l8l[l9|. Sup- 
pose that an arbitrary two-dimensional quantum state is 
available near a particular lattice site, and the task is 
to quantum teleport [3 the state to a neighboring lat- 
tice site, say the one on the other rail. The fidelity of 
teleportation for such an exercise is given by 
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However, the same protocol that tcleports the quantum 
state to the NN site on the other rail, will also teleport 
the state, with a different teleportation fidelity, to two 
other NN sites on the same rail pJl] , with the new fidelity 
being 
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Numerical simulations show that although there are 
two types of bipartite entanglement, which show com- 
plementary behavior with increasing system size, the re- 
gional entanglement, Pavg, decreases as the size of the 



lattice is increased (see Fig. 4). 
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FIG. 4: (Color online.) Regional entanglement in RVB liquids 
on ladders of different sizes. 



To sum up, regional entanglement of RVB liquid 
on ladders is a quantity that mirrors the behavior of 
nearest-neighbor bipartite entanglement of isotropic lat- 
tice RVBs. However, the internal picture of bipartite 
entanglement in RVBs is far richer and depends on the 
geometry of the lattice. We have also observed that in 
the case of RVB ladders without periodic boundary con- 
ditions, the complementary behavior between bipartite 
entanglement of rails and steps remains the same as in 
the case of ladders with periodic boundary conditions. 



B. Analytical estimates of bipartite entanglement 

1. Upper Bounds from Monogamy of Entanglement 

A first estimate on the bipartite entanglement can be 
obtained by using monogamy of entanglement (20j . Let 
us consider an arbitrary site iai say on sublattice A. If 
the lattice is finite, and if periodic boundary conditions 
are not assumed, we assume that ia is not on a bound- 
ary step (see Fig. 1). Each such site is surrounded by 
three NN sites, belonging to the sublattice B. Rotational 
invariance of the state j?/;) ensures that all the three NN 
bipartite states are in Werner states, of which two have 
Werner parameter jOj., and one has Pa- The monogamy 
of entanglement [20] demands that 2r(/5j.) -I- t(ps) < 
'''i:rest( I ■'/'))• t{p) is the "tangle" of the bipartite state 
p, and is defined as max(0,Ai — A2 — A3 — A4), where 
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Ai's are the square roots of the eigenvalues, in decreas- 
ing order, of pp. Here, p = {oy ® ay)p*{ay ® Oy), where 
complex conjugation is with respect to the (Jz eigenbasis 
[13, [2^ ■ And T]^:rest(IV')) is the tangle of the state in 
any bipartition of one site to rest of the sites, and is unity. 
Also T{p{pr/s)) = (Spr/s — 1)^/4, SO that the monogamy 
inequality reads 

(3p, - 1)2/2 + {3ps - 1)V4 < 1, (5) 

The bound on and Ps obtained in the form of the above 
inequality is depicted in Fig. 5. The 2D projection of 
Fig. 5 clearly shows that the allowed value for ps can 
go up to 1 while the same for p^. is < 0.8. Therefore, 
the complementary behavior between the entanglements 
in the rails and steps can also be anticipated from the 
monogamy inequality. 



for all possible 0. Note here that if ^ 0, it would im- 
ply Pr ^ and < 1. Moreover, Pr < will imply 
no bipartite entanglement in the rails while there can be 
some ~ in principle, maximal - entanglement in the steps. 
Numerical evidence (in Figs. 2 and 3) also suggest that 
this is indeed the trend of Pr and Ps- In case of a multi- 
partite state, maximal bipartite entanglement in any of 
its two-party reduced density matrices leads to no gen- 
uine multipartite entanglement. The question remains 
whether 9 will be zero for large system-size. 

To obtain estimates of 9, we insert the values of and 
Ps, that have been obtained from the numerical simula- 
tions (for different N) in Sec. IIIA, in the inequalities 
in ([6|). We then solve the above inequalities to obtain 
9 G for the first inequality for a fixed N, and similarly 
9 e S2 ■ We now consider the allowed 9 lying in the inter- 
section S^nS^, and plot 6'n,ax = max{6l : 61 e 5f n 5^} 
with respect to N. We find that ^max is decreasing with 




Pt < ^{sin^9 + V2sm29); 
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FIG. 5: (Color online.) Estimate of Werner parameters from 
the monogamy of entanglement. The allowed combinations of 
Pr and Ps are those for which the surface (3pr — 1)^/2-1- (Sps — 
1)^/4 — 1 is negative, and can be read off from the projection 
of the surface. 



2. Bounds from Asymmetric Cloning 

Consider again the task of sending an arbitrary auxil- 
iary qubit from any site iA to its three neighboring sites 
[la |. The teleportation fidelity of the output state at a 
site on the same rail (as i^) will be {pr + l)/2, while the 
fidelity will be (p^ -f l)/2 of the output state on the other 
rail .19.] . Quantum mechanics implies that the fidelities 
of the output states in the whole protocol cannot exceed 
the fidelities of three approximate clones of the initial 
state in the optimal 1 — )■ (1 -I- 2) asymmetric quantum 
cloning machine [2l|, [S^). Therefore, we obtain 
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FIG. 6: (Color online.) Decrease of ^max for increasing 
N. The dots represent Smax obtained from Eq. ([6]) cor- 
responding to A'^ = 6, 8, 10, 12. Pink and green lines are 
respectively linear (0.747664 — 0.0185155a;) and quadratic 
(0.671077 - 0.0010471a;^) fits of the points. The correspond- 



ing mean square errors are 1.22 x 10 
1.06 X 10~^ for the quadratic. 



for the linear fit and 



P, < 1 



(6) 



the increase of the size of the lattice, as shown in Fig. 6. 
The inequalities in Eq. ^ also shows the complemen- 
tarity between the entanglement of rails and steps (see 
Fig. 7). 

This gives us further evidence that there is negligible 
or no multipartite entanglement for large lattice size, in 
contrast to the case of isotropic lattices [6] . We will con- 
cretize this evidence by calculating a measure of genuine 
multipartite entanglement for RVB liquids on ladders, in 
the next section. 



IV. NEGLIGIBLE GENUINE MULTIPARTITE 
ENTANGLEMENT 

The numerical and analytical studies on nearest- 
neighbor bipartite entanglement in the preceding section 
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FIG. 7: (Color online.) Complementary behavior of the en- 
tanglement on the rail with that on the step. The bound on 
Pr and Ps, as obtained from the Eq. for all possible val- 
ues of 9, is depicted by the ellipse in the figure. The allowed 
combinations of pr and Ps fall inside the ellipse. The (red) 
open curve is the quadratic fit (-0.858a;^ -I- 0.241a; -I- 0.67) of 
the (red) circles, which are in turn the values found by ex- 
act calculations for N = 6,8, 10, 12. All calculated points are 
within the bounding ellipse. It is evident that as ps attains 
the maximum value 1, the value of pr goes to 0. 



side of the bipartition. This can be explained by the com- 
plementary behavior of bipartite entanglement in steps 
and rails as discussed earlier. Therefore, the trend of 
multipartite as well as bipartite entanglements indicate 
that for large N, only bipartite entanglement in the steps 
will remain, while bipartite entanglement of the rails as 
well as multipartite entanglement of the whole RVB liq- 
uid will disappear. 



already suggest that the RVB liquid on a ladder has neg- 
ligible or no genuine multiparty entanglement. This is 
due to the evidence presented that the entanglement of 
the states on the steps are near-maximal or maximal. A 
maximally entangled state in d(E) d must be pure [23^ (cf. 
|24]), and hence cannot have any correlation, classical or 
quantum, with any other quantum system (cf. [25j|). 

To quantify the amount of genuine multiparty entan- 
glement present in the RVB liquid on ladders of different 
sizes, we consider a recently introduced genuine multi- 
partite entanglement measure called generalized geomet- 
ric measure (GGM) [H [13. The GGM of an A^-party 
pure quantum state |(/>jv) is defined as 

£:(|0Ar)) = l-AL.(|0iv)), (7) 

where Amaxd^Ar)) = niax Kxl^jv) I, with the maximiza- 
tion being over all pure states |x) that are not genuinely 
Af-party entangled. It was shown in Ref. [14] that 

£{\(f>N}) = l-ma.x{X\.B\A\JB = {1,2, . . . ,N},AnB = 0}. 

. . (8) 

where Xa-.b is the maximal Schmidt coefficients in the 
A : B bipartite split of |(^jv) ■ 

In Fig. 8, we find that indeed the genuine multipartite 
entanglement measure £ for the RVB liquid on a ladder, 
decreases with increasing A^. We also notice that when 
performing the maximization in Eq. ([5]) for obtaining 
the GGM, the maximum Schmidt coefiicient is obtained 
when the maximum number of steps is included on one 



FIG. 8: (Color online.) Genuine multipartite entanglement 
measure of RVB liquid on ladders. The figure clearly shows 
that the genuine multipartite entanglement decreases with the 
increase of system size A^. 



V. CONCLUSION 

We have considered the resonating valence bond liq- 
uid on a ladder with periodic boundary conditions. We 
have found two different kinds of bipartite entanglement: 
While the bipartite entanglement on the steps is increas- 
ing, that of the rails is decreasing, with the increase of 
the size of the lattice. Both numerical and analytical 
bounds support this thesis. Moreover, genuine multipar- 
tite entanglement of the ladder decreases with increasing 
system size. This is in sharp contrast with the situation 
in isotropic lattices, where same state (RVB liquid) has 
negligible bipartite entanglement, but substantial mul- 
tipartite entanglement. This shows that geometry can 
play an important role in determining the entanglement 
properties of multiparty quantum states. 



Acknowledgments 

We thank Ujjwal Sen for helpful discussions. 



6 



[1] L. Pauling, Proc. R. Soc. London, Ser. A 196, 343 (1949); 
L. Pauling, The Nature of Chemical Bonds (Cornell Uni- 
versity Press, Ithaca, NY, 1960). 

[2] P. Fazekas and P.W. Anderson, Philos. Mag. 30, 23 
(1974). 

[3] P.W. Anderson, Science, 235, 1196 (1987). 

[4] G. Baskaran, Indian J. Phys. 89, 583 (2006). 

[5] A. Y. Kitaev, Ann. Phys. (Leipzig) 303, 2 (2003). 

[6] A. Chandran, D. Kaszlikowski, A. Sen(De), U. Sen, and 

V. Vedral, Phys. Rev. Lett. 99, 170502 (2007). 
[7] E. Dagatto, J. Riera and D. Scalapino, Phys. Rev. B 45, 

5744 (1992). 

[8] E. Dagatto and T.M. Rice, Science 271, 618 (1996). 
[9] S. R. White, R.M. Noack, D. J. Scalapino, Phys. Rev. 

Lett 73, 886 (1994). 
[10] G. Sierra and M.A. Martin-Delgado, Phys. Rev. B 56, 

8774 (1997). 

[11] A. Tribedi and I. Bose, Phys. Rev. A 79, 012331 (2009; 

LA. Kovacs and F. Igloi, Phys. Rev. B 80, 214416 (2009), 

and references therein. 
[12] Y. Li, T. Shi, B. Chen, Z. Song, and C.-P. Sun, Phys. 

Rev. A 71, 022301 (2005). 
[13] A. Sen(De) and U. Sen, Phys. Rev. A 81, 012308 (2010) 



[14] A. Sen(De) and U. Sen, larXiv:1002. 1253 [quant-ph] 
[15] S. Liang, B. Doucot, and P.W. Anderson, Phys. Rev. 

Lett. 61, 365 (1988). 
[16] R.F. Werner, Phys. Rev. A 40, 4277 (1989). 
[17] W.K. Wootters, Phys. Rev. Lett. 80, 2245 (1998). 
[18] C.H. Bennett, G. Brassard, C. Crepeau, R. Jozsa, A. 

Peres, and W.K. Wootters Phys. Rev. Lett. 70, 1895 

(1993). 

[19] P. Horodecki, M. Horodecki, and R. Horodecki, Phys. 

Rev. A, 60, 1888 (1999). 
[20] V. Coffman, J. Kundu, and W.K. Wootters, Phys. Rev.A 

61, 052306 (2000); T.J. Osborne and F. Verstraete, Phys. 

Rev. Lett. 96, 220503 (2006). 
[21] M. Murao, D. Jonathan, M.B. Plenio, and V. Vedral, 

Phys. Rev. A 59, 156 (1999). 
[22] S. Ibhsdir, A. Acm, N.J. Cerf, R. Filip, J. Fiurasek, and 

N. Gisin, Phys. Rev. A 72, 042328 (2005). 
[23] Private communication with M. Horodecki, U. Sen, and 

W.K. Wootters, of one of the authors. 
[24] R. Garisto and L. Hardy, Phys. Rev. A 60, 827 (1999). 
[25] A. Ekert, Phys. Rev. Lett. 67, 661 (1991). 



